Relativistic tunneling and accelerated transmission 
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Abstract 

We obtain the solutions for the tunneling zone of a one-dimensional electrostatic potential in the 
relativistic (Dirac to Klein-Gordon) wave equation regime when the incoming wave packet exhibits 
the possibility of being almost totally transmitted through the potential barrier. The conditions 
for the occurrence of accelerated and, eventually, superluminal tunneling transmission probabilities 
are all quantified and the problematic superluminal interpretation originated from the study based 
on non-relativistic dynamics of tunneling is reevaluated. The treatment of the problem suggests 
revealing insights into condensed-matter experiments using electrostatic barriers in single- and 
bi-layer graphene, for which the accelerated tunneling effect deserves a more careful investigation. 
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Finding a definitive interpretation for the nature of swperluminal barrier tunneling has 
brought up a fruitful discussion in the literature since pulses of light and 

microwaves appear to tunnel through a barrier at speeds faster than a reference pulse moves 
through a vacuum (], 3, Isi 9]. Tunneling occurs when a wave impinges on a thin barrier 
of opaque material and some small amount of the wave leaks through to the other side. 
The superluminal experiments that promoted the controversial discussions were performed 
with a lattice of layers of transparent and opaque materials arranged so that waves of 
some frequencies are reflected (through destructive interference) but other frequencies pass 
through the lattices in a kind of filter effect correlated to the Hartman effect 101 ]. 

In all cases described by the non-relativistic (Schroedinger) dynamics Ejj], the pulse (wave 
packet) that emerges from the tunneling process is greatly attenuated and front-loaded due 
to the filter effect (only the leading edge of the incident wave packet survives the tunneling 
process without being severally attenuated to the point that it cannot be detected). If one 
measures the speed by the peak of the pulse, it looks faster than the incident wave packet. 
Since the transmission probability depends analytically on the momentum component k 
(T = T(k)), the initial (incident wave) momentum distribution can be completely distorted 
by the presence of the barrier of potential. As there is no sharp beginning to a pulse, we 
cannot declare the instant of its arrival at a certain point. Thus the computation of the 
tunneling time becomes fundamentally meaningless. We could only watch the rising edge of 
the pulse and try to recognize what is arriving 

By employing a tunneling dynamics described by a relativistic wave equation we can 
reevaluate the most part of these difficulties. In this scenario, the existence of formal analo- 
gies between the barrier tunneling of the pulses of light and the tunneling transmission of 
relativistic particles allows for a close correspondence between quantum relativistic motions 
described by the Klein-Gordon equation and electromagnetic wave propagations in the pres- 
ence of dissipation [11]. Here we demonstrate with complete mathematical accuracy that, in 
some limiting cases of the relativistic (Klein) tunneling phenomena where the relativistic ki- 
netic energy is approximately equal to the potential energy of the barrier, and mcL/h << 1, 
particles with mass m can pass through a potential barrier V of width L with transmis- 
sion probability T approximately equal to one. Since T ~ 1, the analytical conditions for 
the stationary phase principle applicability which determines the tunneling (phase) time 
for the transmitted wave packets are totally recovered. Differently from the previous (non- 
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relativistic) tunneling analysis, the original momentum distribution is kept undistorted and 
there is no filter effect. The tunneling time is then computed for a completely undistorted 
transmitted wave packet, which legitimizes any eventual accelerated transmission. 

Some authors consider difficult and perhaps confusing the treatment of all interactions 
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ane waves or wave packets with a barrier potential using a relativistic wave equation 
. This is because the physical content depends upon the relation between the 
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barrier height V and the mass m of the incoming (particle) wave, beside of its total energy 
E. In the first attempt to evaluate this problem, Klein [16j] considered the reflection and 
transmission of electrons of incidence energy E on the potential step V(x) = Q(x)V in the 
2-dimensional time- i ndep endent Dirac equation which can be represented in terms of the 
usual Pauli matrices [17j by [31] 



[a 3 a l di - (E - Q(x 1 )V ) - a 3 m] 4>(k,x 1 ,x 2 ) = 0, (from this point c = h = 1), (1) 

which corresponds to the reduced representation of the usual Pauli-Dirac gamma matrix 
representation (i = 1, 2). The physical essence of such a theoretical configuration lies in 
the prediction that fermions can pass through large repulsive potentials without exponential 



damping. It corresponds to the so called Klein tunneling phenomenon 13j which follows 



accompanied by the production of a particle-antiparticle pair inside the potential barrier. 
It is different from the usual tunneling effect since it lies in the energy zone of the Klein 
paradox [if], l^ - Taking the quadratic form of the above equation reduced to one-dimension 
for a generic scalar potential V(x), we obtain the analogous Klein-Gordon equation, 

(ido - V{x)f 4>{k, x) = (E- V{x)f cj){k, x) = {-d 2 x + m 2 ) cf)(k, x), (2) 

which, from the mathematical point of view, due to the second-order spatial derivatives, 
has similar boundary conditions to those ones of the Schroedinger equation and leads to 
stationary wave solutions characterized by a relativistically modified dispersion relation. 

By depicting three potential regions by means of a rectangular potential barrier V(x), 
V(x) = V if < x < L, and V(x) = if x < and x > L, we observe that the incident 
energy can be divided into three zones. Differently from the energy configuration relative 
to the non-relativistic (Schroedinger) dynamics, the above barrier energy zone, E > V + m, 
involves diffusion phenomena of oscillatory waves (particles). In the so called Klein zone 



m 



16j, E < V — m, we find oscillatory solutions (particles and antiparticles) in the barrier 



region. In this case, antiparticles see an opposite electrostatic potential to that seen by 



the 
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aarticles and hence they will "see" a potential well where the particles "see" a barrier 
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18, 19]. The tunneling zone, V — m < E < V + m, for which only evanescent waves exist 



in the barrier region, is that of interest in this work. 
By evaluating the problem for this tunneling (evanescent) zone assuming that <j)(k, x) are 
stationary wave solutions of the Eq. (j2J), when the peak of an incident (positive energy) wave 
packet reach the barrier x = at t = 0, we can write 

4>i(k,x) = exp [i k x] + R(k, L) exp [— i k x] x<0, 
(j) 2 {k,x) = a(k) exp [— p(k) x] + /3(k) exp [p(k) x] < x < L, (3) 
(j) 3 (k,x) = T(k, L) exp [i k{x — L)] x > L, 

where the novel dispersion relations: k 2 = E 2 — m 2 and p(k) 2 = m 2 — (E — V ) 2 , are modified 
with respect to the usual non-relativistic ones. 

In order to proceed with a phenomenological analysis which allows us to establish a corre- 
spondence with the non-relativistic (NR) solutions, it is convenient to define the kinematic 



<p{k,x) 



variables in terms of the following parameters: w = V2mV^, v = V /m = w 2 /2m 2 } and 
n 2 (k) = k 2 /w 2 = Enr/V . The parameter w corresponds to the same normalization param- 
eter of the usual NR analysis where k 2 = 2mE^R- The previously quoted relation between 
the potential energy V and the mass m of the incident particle is given by the parameter 
v. Finally, n 2 (k) represents the dependence on the energy for all the results that are being 
considered here. After simple mathematical manipulations, it is easy to demonstrate that 
the tunneling zone for the above Klein-Gordon equation (J2J) is comprised by the interval 
(n 2 (k) — v/2) 2 < 1 which made n 2 (k) assume larger values {n 2 (k) » 1), in opposition to 
the NR case where the tunneling energy zone is constrained by < n 2 (k) < 1). We shall 
observe that such a peculiarity has a suitable relation with the possibility of superluminal 
transmission through the barrier. The limits for NR energies (k 2 « m 2 and V « m) 
given by vn 2 (k) « 1 and v/n 2 (k) << 1 reproduce the Schroedinger equation results for the 
transmission coefficient and for the corresponding traversal time. 

With regard to the standard one-way direction wave packet tunneling, for the set of sta- 
tionary wave solutions given by Eq. it is well-known 
T(n,L) = \T(n,L)\ exp [i(f(n, L)\ is written in terms of 

\T(n,L)\ = \l + - -L smh 2 [p(n)wL]} . ( i) 



231 ] that the transmitted amplitude 



where we have suppressed from the notation the dependence on k, and 

if(n, L) = arctan < —— tanh \p(n) wL] > , (5) 

I 2np{n) ' J 

for which we have made explicit the dependence on the barrier length L (parameter wL), 



and we have rewritten p(k) = wp(n), with p(n) 2 = y/T+ 2n 2 v — (n 2 — v/2). 

The additional phase <£>(n, L) that goes with the transmitted wave is utilized for calcu- 
lating the transit time t v of a transmitted wave packet when its peak emerges at x — L, 

dk dn(k) d(p(n, L) (L) 1 dip(n,L) 
t(p ~ dE(k) ~dk ~ ~w(L) ~d7i ' ' 

evaluated at k = k (the maximum of a generic symmetrical momentum distribution g(k — k ) 

that composes the incident wave packet). By introducing the classical traversal time defined 

as = L(dk/dE(k)) = L/v, we can obtain the normalized phase time, 

(7) 



V _ f(n,L) 



f(n, L) = 8n 2 [(2 + 8n 2 v + v 2 ) - (in 2 + 3v) y/l + 2n 2 



+ 4 [(4 + An 2 v + v 2 ) Vl + 2n 2 i; - 2v (2 + 3n 2 v) 



s'mh(p(n)wL) cosh(p(n)wL) 



p(n)wL 

g(n,L) = 16n 2 [2 (1 + 2n 2 v) - Vl + 2n 2 v (2n 2 + v)] 

+ 2 [(4 + 8n 2 v + v 2 ) VI + 2n 2 v - Av (1 + 2n 2 v)] smh(p(n)wL) 2 . 

We compare the theoretical results for the tunneling phase times in correspondence with 
their respective transmission probabilities in the Fig.([T]) for different propagation regimes 
(v = 0(NR), 1, 2, 5, 10)). It is important to emphasize that the tunneling region is com- 
prised by the interval (n 2 — v/2) 2 < 1, n 2 > 0. 

We can notice the possibility of accelerated (t v < t^), and eventually superluminal 
(negative tunneling delays, t 9 < 0) transmissions without recurring to the usual analysis 
of the opaque limit (p(n)wL — > oo) which leads to the Hartman effect 10J. In the NR dy- 
namics (Schroedinger equation solutions), the opaque limit and its consequent superluminal 
interpretation (Hartman effect) were extensively discussed in the literature. Superluminal 
group velocities in connection with quantum (and classical) tunneling were predicted even 
on the basis of tunneling time definitions more general than the simple Wigner's phase-time 
24] (Olkhovsky et al, for instance, discuss a simple way of understanding the problem Q]). 
Discussions on relativistic causality in addition to several analytical limitations have ruined 
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Tunneling 
Klein 
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Normalized Energy n (k) = k 2 /2mV 

FIG. 1: Tunneling TRANSMISSION probabilities and the corresponding tunneling PHASE times 
for the dynamics of the relativistic wave equation. We have classified the energy zones by the 
line thickness: the thick line corresponds to the tunneling energy zone, the intermediate line 
corresponds to the Klein zone and the thin line corresponds to the above barrier energy zone. In 
fact the tunneling region is comprised by the interval intersection (n 2 — v/2) 2 < 1 and n 2 > 0. Here 
we have adopted the illustratively convenient value of wL = 2ir. We have set v = 0, 1, 2, 5, 10 and 
we have constrained our analysis to n 2 (k) > since we have assumed V > 0. It is convenient to 
observe that the NR regime can be parameterized by the limit where v — > 0. 
the mo st p art of possibilities of superluminal interpretation of the NR tunneling phenom- 
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231 ] . In a causal manner, the arguments consist in explaining the superluminal 



phenomena during tunneling as simply due to a reshaping of the pulse, with attenuation, as 
already attempted (at the classical limit) [251 ]. i. e. the later parts of an incoming pulse are 
preferentially attenuated, in such a way that the outcoming peak appears shifted towards 
earlier times even if it is nothing but a portion of the incident pulse forward tail [?], 0]. 

We do not intend to extrapolate to the delicate question of whether superluminal group- 
velocities can sometimes imply superluminal signalling, a controversial subject which has 
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been extensively explored in the literature ([5| and references therein). Otherwise, the 
phase time calculation based on the relativistic dynamics introduced here offers distinct 
theoretical possibilities of superlunminal transmission in a novel scenario, for the limit case 
where p(n) wL tends to (with L ^ 0), in opposition to the opaque limit where p(n) wL tends 
to oo. Let us then separately expand the numerator f(n,L) and the denominator g{n,L) of 
the Eq. (j7|) in a power series of p(n) wL (p(n) — ► 0) in order to observe that in the lower 
(upper) limit of the tunneling energy zone, where n 2 tends to v/2 + (— )1, the numerical 
coefficient of the zero order term in p(n) wL amazingly vanish in the numerator as well as in 
the denominator! Since the coefficient of the linear term also is null, just the coefficient of 
the second order terms plays a relevant role in both series expansions. After expanding the 
Eq. ([7j), such a step-by-step mathematical exercise leads to 

t v 4 [(4 + An 2 v + v 2 ) VI + 2n*v - 2v (2 + 3n\>)] 

— — = — f , f + C/(p(n) 8 

r (fc) 3 [(4 + 8n 2 v + u 2 ) VI + 2n 2 t; - 4u (1 + 2n 2 ^)] 

for small values of p(n). At the same time, since \im n 2^ v / 2 ^i p(n) = 0, the tunneling trans- 
mission probability can be approximated by 



lim \T(n,L) 



1+ (wL) 2 



2u=f4 



w » i 

2 



[1 + (mL) 2 H , (9) 



from which we recover the high probability of complete tunneling transmission when mL < < 
1. Finally, for the corresponding values of the phase times evaluated in (JSJ), we obtain, 



lim — p — = , n 2 ^u/2=Fl, n 2 , v > 0, (10) 

t/m 31±2n 2 ' v 7 

that does not depend on mL, and we notice that its asymptotic (ultrarelativistic) limit always 

converges to 0. In particular, in the lower limit of the tunneling energy zone, n 2 — > v/2 — 1, 

it is always negative. Since the result of Eq. ffTOl) is exact, and we have accurately introduced 

the possibility of obtaining total transmission (transparent barrier), our result ratifies the 

possibility of accelerated transmission (positive time values), and consequently superluminal 

tunneling (negative time values), for relativistic particles when mL is sufficiently smaller than 

1 (=^» T « 1). Keeping in mind that the barrier height has to be chosen such that one remains 

in the tunneling regime, it is notorious that the transmission probability depends only weakly 

on the barrier height, approaching the perfect transparency for very high barriers, in stark 

contrast to the conventional, non-relativistic tunneling where T(n,L) exponentially decays 



with the increasing V . Since this results can be analytically extended to the Klein paradox 
energy zone, such a relativistic effect is usually attributed to a sufficient strong potential 
that, being repulsive for electrons, is attractive for positrons and results in unstable positron 



states inside the barrier, which align the energy with the electron continuum outside [26]. 

Obviously, the above results correspond to a theoretical prediction, in certain sense, not 
so far from the experimental realization. By considering the magnitude of the parameter 
mL (mc 2 /[h(c/ L)] in standard units) for an electron with mass ~ 0.5 Mel 7 , and observing 
that in natural units we have 0.2 MeV pm ~ 1, we conclude that it should be necessary a 
potential barrier of width L « lpm to permit the observation of the quoted superluminal 
transmission. By principle, its observation makes the effect relevant only for some exotic 
situations as positron production around super-heavy nuclei (Z ~ 170) [27] or evaporation 



of black holes through generation of particle-antiparticle pairs near the event horizon [28]. 

In the most common sense, the above condition should be naturally expected since we 
are simply assuming that the Compton wavelength (h/(mc)) is much larger than the length 
L of the potential barrier that, in this limit situation, becomes invisible for the tunneling 
particle. The relativistic quantum mechanics establishes that if a wave packet is spread 
out over a distance d » 1/m, the contribution of momenta \p\ ~ m » 1/d is heavily 
suppressed, and the negative energy components of the wave packet solution are negligible; 
the one-particle theory is then consistent. However, if we want to localize the wave packet 
in a region of space (wave packet width d) smaller than or of the same size as the Compton 
wavelenght, that is d < 1/m, the negative energy solutions (positron states) start to play an 
appreciable role. This qualitative arguments report us to the Klein paradox and the creation 
of particle-antiparticles pairs during the scattering process which might create the intrin- 
sic (polarization) mechanisms for accelerated and/or non-causal fermion teletransportation. 
The condition d < L < 1/m (where d < L is not mandatory) imposed over a positive energy 
component of the incident wave packet in the relativistic tunneling configuration excite the 
negative energy modes (antiparticles) in the same way that the movement of electrons in a 
semi-conductor is concatenated with the movement of positively charged holes. 

Turning back to the context of the nanoscopic scale structures, the most challenging pos- 
sibility of observing similar effects occurs for massless (or effective mass) Dirac fermions in 
graphene structures. Even though the linear spectrum of fermions in graphene implies zero 
rest mass, their cyclotron mass approaches to 10~ 2 m e [29], which increases the superluminal 
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tunneling scale to 1 angstrom. In spite of the theoretical focus, the results here obtained 
apply to some configurations which should deserve further attention by experimenters in 
the study of the graphene structures where the dynamics of the electron is described by a 
relativistic-like dynamics. For bi-layer structures, due to the chiral nature of their quasiparti- 
cles, quantum tunneling in these materials becomes highly anisotropic, qualitatively different 
from the case of normal, non-relativistic electrons. In fact, it has been recently speculated 
that, from the experimental point of view, the graphene provides an effective medium for 
mimicking relativistic quantum effects where, for instance, massless Dirac fermions allow 
a close realization of Kleins gedanken experiment whereas massive chiral fermions in bi- 
layer graphene offer an interesting complementary system that elucidates the basic physics 
involved. The point is that, in conventional two-dimensional systems, sufficiently strong 
disorder results in electronic states that are separated by barriers with exponentially small 
transmittance 26|]. In contrast, in single- and bi-layer graphene materials all potential bar- 



riers are relatively transparent (T(n,L) ~ 1): the quasiparticles in graphene exhibit a linear 
dispersion relation E = hkvf that corresponds to the pseudo-ultrarelativistic limit of our 
analysis for pseudo- massless particles traveling with their Fermi velocity Vf. In this case 
there are pronounced transmission resonances where T approaches unity for some particular 
geometric configurations, which does not allow charge carriers to be confined by potential 
barriers that are smooth on atomic scale. Moreover, some authors have demonstrated ex- 
perimentally and theoretically, that the biased graphene bilayer is a tunable semiconductor 
where the electronic gap can be controlled by the electric field effect reaching values as large 
as 0.3 eV Q. 

To summarize, as previously pointed out, we have considered tunneling by a wave packet 
obeying the Klein-Gordon equation, in a regime where the Compton wavelength of the 
particle is much larger than the width of the barrier, assumed rectangular. Our standard 
analysis have suggested the possibility of total transmission together with an accelerated 
tunneling. Once it has been observed that the physical essence of the Klein paradox lies in 
the prediction that particles can pass through large repulsive potentials without exponential 
damping [13I . [14J] , the scenario of the Klein-paradox and accelerated tunneling transmission 
associated with relativistic-like phenomena at nanoscopic scale can be tested experimentally 
using graphene devices. In a subsequent study we intend to investigate the appearance of an 
equivalent smaller effective mass value M e // << m due the minimal coupling of the charged 
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particle magnetic momentum with an external magnetic field, in particular, for the massless 
electron propagation in a single-layer graphene which could introduce some novel ingredients 
for quantifying these peculiarities of the relativistic tunneling effect. 
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